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ABSTRACT 


The  effectiveness  of  the  iterative  beaming  technique  for  the  separation  of 
superposed  Rayleigh  wave  signals  was  evaluated.  This  technique  can  be  useful 
for  the  detection  and  discrimination  of  seismic  events  the  signals  of  which 
arrive  simultaneously  at  a seismic  array. 

The  technique  was  applied  to  superposed  recordings  of  Rayleigh  waves  at 
LASA,  It  was  found  that  multipathing  of  Rayleigh  waves  can  cause  enough  leakage 
into  the  beam  of  the  desired  event  to  make  the  separation  impossible  if  the 
amplitude  of  waves  from  the  event  to  be  eliminated  is  more  than  three  times 
that  of  the  event  of  interest. 
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INTRODUCTION 


It  has  been  shown  theoretically  that  the  iterative  beam  processor  is 
equivalent  to  the  maximum  likelihood  processor  (Blandford  et  al. , 1973; 

Shumway,  1972;  Dean  et  al.,  1968).  The  method  has  been  successfully  demon- 
strated on  short  period  body  wave  data  at  TFO  and  was  shown  to  be  superior 
to  simple  beaming.  This  report  evaluates  the  performance  of  the  method 
when  it  is  applied  to  mixed  surface  waves. 

The  iterative  beam  processor  assumes  that  the  directions  of  the  two 
sources  are  known.  By  beaming  the  array  at  the  first  epicenter  one  obtains 
a first  signal  estimate  for  the  first  event  (which  is  contaminated  by  the 
second  event).  The  time  shifted  first  signal  estimate  is  then  subtracted  from 
all  sensor  traces  and  the  resulting  traces  are  beamed  on  the  second  event  to 
obtain  first  estimate  of  the  second  signal.  This  estimate  is  also  subtracted 
from  the  original  traces  and  the  result  is  rebeamed  on  the  first  event  to 
obtain  a second  estimate  for  the  first  signal.  Repetition  of  this  procedure 
should  theoretically  yield  maximum  likelihood  estimates  for  the  two  signals 
(Blandford  et  al.,  1973).  The  computation  involves  only  time  shifting  and 
summing  instead  of  the  complicated  matrix  operations  and  convolutions  needed 
in  computing  maximum  likelihood  filters  tue  usual  way.  If  the  directions  of 
the  two  (or  more)  signal  arrivals  are  not  known  the  iterative  beaming  method 

cannot  be  used. 

Although  the  method  can  be  applied  to  long  period  body  waves,  long  period 
surface  waves  are  more  important  in  discrimination  studies  and  we  shall  apply 
the  method  to  superpositions  of  Rayleigh  waves.  The  method  was  shown  to  be 
theoretically  effective  for  superposed  plane  waves,  but  it  has  to  be  tested 


TR-73-7,  Teledyne  Geotech,  Alexandria,  Virginia. 
Report  it 207,  Teledyne  Geotech,  Alexandria,  Virginia. 


-7- 


tual  surface  wave  data  since  surface  waves  are  commonly  refracted  lateral!,, 
resulting  In  multipathing  and  deviations  of  the  directions  of  preach  of  the 
signal  from  the  great  circle  path. 

The  usefulness  of  the  method  will  be  limited  by  these  properties  oi 

, v.  nf  these  limitations  should  be  determined 

surface  waves,  and  the  extent  of  these  1 iteratlve  bea.s  ln  the 

Pinoi  rically  The  amount  of  energy  leaking 

empirically.  Th  direction  impairs  the  detectability  of 

directions  other  than  the  epicent  since  the 

surface  waves  from  weak  events  in  the  presence  of  a large 

surface  wave  magnitude  is  an  important  discriminant  these  limrtatrons  must  be 

The  method  is  also  suitable,  as  shown  below,  to  quickly  separu 
known.  The  metnoa  offprtive  research 

multipath  arrivals  in  the  time  domain  and  can  be  used  as  an  effective 

tool  for  the  understanding  of  multipathing  in  specific  cases. 


DATA 


A Get  of  events  which  was  previously  used  in  a study  of  frequency  wave- 
number spectra  (Lambert  and  Der,  1974)  was  utilized  in  this  report.  The 
epicentral  data  of  these  events  are  given  in  Table  I.  Vertical  component 
long  period  seismometer  recordings  at  LASA  were  used  in  this  study.  Both  the 
full  array  and  a part  of  the  full  array  consisting  of  the  center  element  Aq 
and  the  B and  C rings  was  used  in  the  analysis.  Two  elements  of  C rings 
(Cl  and  C2)  were  omitted  due  to  data  problems.  In  addition,  the  element  D was 
omitted  for  the  Argentina  event  for  the  same  reason.  In  addition  to  individual 
event  recordings,  linear  combinations  of  seismic  traces  from  several  events 
and  noise  were  also  analyzed.  The  same  combinations  were  used  in  a previous 
study  (Lambert  and  be*. , 19  74). 


Lambert,  J.  and  Der,  Z.  A.,  1974.  Comparison  of  two  segment  maximum  likeli- 
hood (TSML)  frequency  wavenumber  spectra  with  the  fast  beamed  frequency 
wavenumber  spectra  (FKPLOT) , Seismic  Data  Analysis  Center  TR-74-6, 
Teledyne  Geotech,  Alexandria,  Virginia. 


DATA  PROCESSING 


Recordings  of  individual  events  were  processed  by  the  iterative  beaming 
technique  by  assuming  various  directions  for  a second  (absent)  signal  and 
examining  the  leakage  of  the  energy  from  the  event  into  the  iterative  beam. 
Iteration  is  stopped  if  the  rms  difference  between  successive  iterations  was 
less  than  2%  or  if  the  number  of  iterations  exceeded  six.  This  assumes  that 
the  above  criteria  are  sufficient  to  assure  that  the  theoretical  limit  of 
improvement  was  sufficiently  approximated.  The  ratio  of  the  rms  signal  relative 
to  the  rms  iterative  leakage  beam  and  the  simple  leakage  beam  was  plotted  on 
polar  diagrams  in  dB  units  in  order  to  compare  the  relative  effectiveness  of 
the  two  processes.  The  azimuths  of  the  second  signal  were  assumed  to  be  at 
30°  intervals  from  the  azimuth  of  the  first  (single  event)  signal.  Figures  1 
through  6 clearly  show  the  superiority  of  the  iterative  technique  over  the 
simple  beam.  These  diagrams  are  comparable  to  th i f-k  spectra  using  FKPLOT 
and  the  maximum  likelihood  spectra  (Lambert  and  Der,  1974).  It  must  be 
pointed  out  however  that  the  iterative  beaming  is  a wide  band  process  while 
the  f-k  spectra  were  computed  for  a single  frequency.  Nevertheless,  both 
kinds  of  results  demonstrate  the  higher  resolution  property  of  the  maximum 
likelihood  process.  In  discussing  the  results  we  shall  occasionally  refer  to 
the  f-k  spectra  given  in  the  previous  report  (Lambert  and  Der,  1974)  which 
analyzes  the  same  events. 

Combined  events  were  processed  to  produce  signal  estimates  for  both 
events  which  can  be  visually  compared  to  the  original  signals. 
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RESULTS 


In  the  following  we  shall  describe  the  polar  plots  of  energy  rejection 
of  simple  beams  and  iterative  beams  centered  at  directions  other  than  that  of 
the  epicenter  for  various  events.  In  these  plots  the  direction  of  epicenter 
is  marked  by  an  arrow.  The  rejection  in  this  direction  is  naturally  zero. 

In  the  other  directions  spaced  at  30°  intervals  the  ratio  of  the  rms  amplitudes 
of  the  event  (averaged  over  the  various  sensors)  and  that  of  the  beam  and 
iterative  beam  are  plotted  in  decibel  units. 

Figure  1 shows  these  plots  for  the  full  LASA  array  obtained  for  the  Baja 
California  event.  This  even,  has  a short  surface  wave  train  not  exceeding  in 
duration  the  window  length  (256  seconds)  used.  Figure  7d  shows  the  locations 
for  three  overlapping  time  windows.  The  middle  one  includes  the  full  surface 
wave  train,  the  first  one,  the  first  half  and  the  third  one,  the  second  half 
of  it.  Beams  have  been  computed  for  three  velocities  3.34,  3.56,  and  3.79 
km/sec.  For  the  first  window  the  best  rejection  is  obtained  at  the  two  higher 
velocities,  for  the  second  and  third  at  the  two  lower  velocities.  This 
undoubtedly  reflects  the  dispersion  of  surface  waves,  the  longer  period  waves 
at  the  beginning  of  the  wavetrain  having  a higher  phase  velocity.  The  velocity 
3.56  km/sec  seems  to  be  a good  overall  beam  velocity  in  this  case  and  it  was 
used  exclusively  for  time  domain  analysis  of  signals  presented  in  the  latter 
part  of  this  document. 

Figure  2 shows  the  results  for  the  reduced  array  containing  only  the 
A,  B,  and  C rings  (as  described  above).  The  dependence  of  the  results  on  the 
beam  velocities  and  the  position  of  the  window  is  about  the  same  as  in  the 
previous  case.  Note  the  additional  rejection  of  off  beam  energy  by  the 
iterative  process,  when  compared  to  the  simple  beam.  An  explanation  for  this 
is  that  the  first  beam  due  to  the  small  array  is  a poor  approximation  to  the 
signal  wavefo  -m  which  improves  greatly  as  the  iteration  proceeds.  Moreover, 
lateral  inhomogeneities  effect  the  wave  shapes  less  within  the  smaller  array 
and  the  wave  appears  to  be  more  coherent.  Since  iterative  beaming  (or  its 
equivalent  the  maximum  likelihood  filtering)  is  critically  dependent  on  the 
assumed  signal  model,  which  is  a coherent  plane  wave,  this  process  will  be 
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markedly  more  efficient  using  the  smaller  array.  The  improvement  of  the 
iterative  technique  over  the  simple  bean  ranges  between  10-15  dB  in  this  case. 

Figure  3 shows  a similar  set  of  plots  for  the  /rgentine  event  for  the 
full  array.  The  location  of  the  data  windows  for  these  diagrams  is  shown  in 
Figure  7e.  These  are  nonoverlapping  windows  sampling  various  parts  of  the 
Rayleigh  wave  train.  The  performance  of  the  iterative  beam  relative  to  the 
beam  is  poor  in  the  average.  Tiie  polar  plots  for  the  reduced  array  show 
visible  improvement  relative  to  the  full  array  beams  (Figure  4).  The 
beginning  of  the  wave  train  yields  better  results  than  the  succeeding  data 
windows.  We  assume  that  multipathing  is  the  cause  of  these  difficulties. 

The  visual  similarity  of  traces  is  poor  and  the  f-k  spectra  obtained  for  this 
event  in  a previous  study  (Lambert  and  Der,  1974)  indicates  multipathing.  We 
shall  show  later  that  a multipath  arrival  in  the  latter  part  of  the  wave  train 
can  actually  be  isolated  by  the  iterative  beaming  technique.  This  opens  up  a 
new  application  of  this  technique  for  studying  multipathing.  The  performance 
of  the  process  changes  less  with  the  beam  velocity  than  in  the  previous 
example,  possibly  because  spectral  components  cf  varying  phase  velocities  are 
mixed  by  multipathing  within  the  same  windows. 

Figures  5 and  6 show  similar  polar  plots  for  the  Guatemala  event.  We 
note  essentially  the  same  features,  slight  dependence  on  velocity  attributable 
to  the  varying  frequency  content  of  the  time  windows,  superior  performance 
of  the  iterative  processor  and  the  reduced  array  outperforming  the  full  array. 
The  location  of  the  time  windows  are  shown  in  Figure  7f.  The  performance 
deteriorates  with  time  indicating  that  the  coda  consists  of  multipathed 
incoherent  arrivals. 

Separation  of  Combined  Events 

Although  the  previous  results  show  how  effective  the  separation  can  be  in 
the  root  mean  square  sense  it  is  desirable  to  see  how  well  the  waveforms  of 
earthquakes  can  be  separated  in  the  time  domain.  If  the  separated  waveforms 
cannot  be  recognized  as  a waveform  from  an  earthquake,  that  is,  dispersed 
surface  wave  trains,  it  would  be  unwise,  for  instance,  to  determine  surface 
wave  magnitudes  from  them  since  they  could  be  multipath  arrivals  or  long 
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period  microseisms.  In  the  following  we  shall  show  the  separation  of  super- 
posed pairs  of  events.  The  seismograms  of  the  individual  events  are  shown  m 
Figures  7a-7f.  The  pairs  of  events  were  superposed  with  roughly  equal 
maximum  amplitudes,  and  decomposed  by  iterative  beaming  into  the  individual 
events.  To  test  for  leakage  of  each  event  into  the  iterative  beam  arrived  at 
the  other  event  the  same  procedure  was  performed  in  the  case  when  only  one 
event  was  present. 


Throughout  the  following  computations  the  two  events  were  adjusted  and 
their  maximum  amplitudes  within  the  time  interval  processed  were  equal.  The 
results  of  the  computations  which  followed  were  plotted  on  a fixed  scale  and 
the  amplitudes  of  the  plots  are  directly  comparaole.  This  was  done  in  order 
to  compare  the  amplitudes  of  the  separated  waveforms  and  the  energy  leakage. 

Figures  8a  and  b show  the  results  of  such  analysis  for  the  superposition 
of  the  Philippine  Islands  event  with  the  New  Britain  event  for  the  full  LASA 
array.  The  traces  on  the  top  of  figure  8a  show  the  results  for  the  Philippine 
event  when  it  was  present  only.  The  leakage  of  energy  into  the  beam  of  the 
New  Britain  event  is  hardly  reduced  beyond  the  first  iteration.  A similar 
procedure  for  the  New  Britain  event  as  shown  on  top  of  figure  8b  shows  visible 
reduction  of  leakage  into  the  beam  of  the  Philippine  event  (traces  1 vs.  2 on 
the  left)  by  about  a factor  of  2 on  the  second  iteration. 


The  middle  traces  in  Figures  8a  and  b are  biased  beams  for  the  individual 
events.  These  traces  are  included  for  comparison  of  waveforms.  Theoretically 
their  amplitudes  are  biased,  less  than  the  actual  amplitudes  of  the  beams  on 
the  individual  events.  The  difference  is,  however,  very  small  only  about  5% 
in  the  worst  case,  the  full  array,  as  calculated  from  the  bias  formula  of 
Blandford,  Cohen,  and  Woods  (1973).  Since  we  need  only  a rough  comparison 
between  the  amplitudes  of  the  beam  and  the  leaked  waveforms  we  can  use  these 
for  comparison.  Unbiased  beams  for  the  same  case  were  not  computed  since 
this  would  have  required  a separate  iteration  sequence. 


As  seen  by  comparing  wavr Forms,  the  ratio  of  the  original  event  amplitudes 
to  that  of  the  leakage  into  the  beam  aimed  at  the  other  event  is  about  4:1 
for  the  Philippine  event  and  7:1  for  the  New  Britain  event.  These  would  be 
roughly  the  limits  imposed  by  multipathing  and  the  array  configuration  in  the 
given  case  in  detecting  small  event  hidden  in  the  tail  of  the  dominant  event. 
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Ambient  noise,  being  of  very  low  amplitude,  does  not  influence  the  results 
in  this  case.  The  bottom  traces  on  Figure  8a  and  b show  the  results  of  the 
decomposition  when  both  events  are  present.  The  successive  iterations  show 
essentially  the  same  waveforms  indicating  that  the  first  beams  of  the  large 
array  achieved  a quite  efficient  separation. 

The  results  of  an  identical  procedure  for  the  reduced  LASA  array  (A,  C, 
and  D rings)  are  shown  in  Figure  9a  and  b for  the  same  two  events.  Because 
of  the  smaller  size  of  the  array  the  iteration  achieves  more  improvement 
relative  to  the  simple  beam.  Nevertheless,  the  leakage  problem  is  more  severe, 
the  amplitudes  (maxima  measured)  of  the  original  amplitude  relative  to  the 
leakage  is  about  2:1  for  the  Philippine  event  and  a somewhat  better  4:1  for 
the  New  Britain  event.  More  iterations  were  also  needed  to  achieve  the  results 
although  not  much  change  can  be  seen  beyond  the  third  unbiased  iteration  for 
any  of  the  events.  The  bottom  traces  again  show  the  separation  results  when 
both  events  are  present.  Since  the  leakage  is  still  fairly  small  not  much 
change  is  seen  in  the  waveforms  of  successive  iterations.  The  azimuthal 
separation  of  the  two  events  in  this  case  was  quite  good,  about  42°,  way 
beyond  the  resolution  capability  of  the  maximum  likelihood  (iterative  beam) 
process  as  seen  in  the  report  by  Lambert  and  Der,  1974.  It  seems  therefore 
that  the  cause  of  "leakage"  is  multipathing,  real  energy  coming  from  the 

] 

undesired  event  in  the  direction  of  the  event  to  be  looked  at.  This  limits 
the  application  of  this  technique  to  superposed  Rayleigh  waves,  and  this 
limitation  will  remain  until  the  multipathing  phenomenon  is  better  understood 
and  is  properly  corrected  for.  Using  a larger  array  reduces  leakage  due  rj 
its  higher  resolution  but  does  not  eliminate  it  completely. 

The  following  figures  show  the  results  of  identical  procedures  applied 
to  other  event  pairs.  Figures  10a  and  b show  the  results  for  the  Philippine 
event  as  combined  with  the  New  Hebrides  event  when  the  full  LASA  array  is 
used.  These  events  are  separated  in  azimuth  by  57°.  Shortly  summarizing 
the  ratio  of  the  original  amplitudes  to  the  leakage  is  about  3:1  for  the 
Philippine  event  and  somewhat  better  than  2:1  for  the  New  Hebrides  event, 
thus  in  spite  the  large  azimuthal  separation  and  the  full  array  used  the 
leakage  is  not  negligible.  As  in  the  previous  case  the  first  beam  is  quite 
effective,  and  no  great  change  in  the  waveforms  can  be  seen  beyond  the  second 


unbiased  iteration.  It  is  a little  hard  to  believe  that  so  much  energy  could 
be  diverted  from  the  main  circle  path,  but  in  the  case  of  New  Hebrides  event 
the  later  arrival  of  the  leaked  energy  seems  to  agree  with  the  notion  that  it 

arrived  over  a longer  path. 

Figures  11a  and  b show  the  results  for  the  reduced  array.  The  ratio  of 
original  to  leaked  wave  amplitudes  is  about  3:1  for  the  Philippine  event  and 
somewhat  larger  than  2:1  for  the  New  Heorides  event,  the  same  as  before.  The 
rest  of  the  remarks  above  also  apply  in  this  case. 

The  next  case  shown  in  Figures  12a  and  b is  the  separation  of  the  Baja 
California  and  Guatemala  events  with  the  full  LASA  utilized.  The  angular 
separation  of  the  azimuths  of  the  epicenters  to  LASA  is  about  43°.  The  ratio 
of  the  original  trace  amplitudes  to  that  of  the  residual  trace  is  about  4:1 
for  the  Baja  event  and  4:1  for  the  Guatemala  event.  The  separation  of  wave- 
forms is  quite  good,  although  some  leakage  of  the  Guatemala  event  into  the 
beam  aimed  at  Baja  is  evident. 

Figures  13a  and  b show  the  same  two  events  processed  using  the  reduced 
LASA  array  only.  The  ratio  of  original  and  leaked  amplitudes  is  about  5:1 
for  the  Baja  event  and  4:1  for  the  Guatemala  event.  Again  more  iterations 
are  needed  to  achieve  convergence.  On  the  separated  superposed  traces  energy 
leaked  from  the  Guatemala  event  is  quite  evident. 

The  results  for  the  event  pair  Baja  California-Argentina  earthquakes 
using  the  full  LASA  array  are  shown  in  Figures  14a  and  b.  The  azimuthal 
separation  of  these  events  is  48°  but  f-k  plots  indicate  that  the  Argentine 
event  is  considerably  multipathed  and  most  of  the  energy  does  not  come  from 
the  actual  azimuth  of  the  source  (Lambert  and  Der,  1974).  Since  the  above 
example  was  calculated  using  the  actual  azimuth  along  the  great  circle  the 
results  are  correspondingly  poor.  The  Argentine  event  leaks  a considerable 
amount  of  energy  into  the  beam  aimed  at  Baja  California,  but  the  beam  aimed 
at  the  actual  azimuth  of  this  event  also  contains  considerable  energy.  The 
separation  of  events  therefore  is  not  successful. 

Figures  15a  and  b show  the  results  for  this  event  pair  using  the  reduced 
arrav,  leading  to  similar  conclusions  as  above. 


-16- 


CONCLUSIONS 


Although  the  iterative  beaming  process  has  been  shown  theoretically  to 
be  very  effective  In  separating  multiple  events  its  application  to  tele- 
seismic  Rayleigh  waves  is  limited  by  the  deviation  of  these  waves  from  the 
assumed  signal  model,  the  plane  wave  with  a constant  phase  velocity.  Wave 
dispersion  which  causes  a gradual  distortion  of  the  plane  waveform  across  the 
array  does  not  seem  to  be  a major  factor  in  reducing  the  effectiveness  of 
this  processor  even  for  large  arrays  (LASA) . The  main  cause  of  the  difficulties 
seems  to  be  multipathing  and  deviation  from  the  great  circle  azimuth  due  to 
lateral  refraction.  The  performance  of  this  processor  depends  on  the  array 
configuration  and  the  location  of  the  events  to  be  separated.  Nevertheless, 
the  above  examples,  although  limited  in  number,  seem  to  indicate  that,  for 
many  cases,  the  separation  cannot  be  relied  on  if  great  circle  a-priori 
azimuths  are  used  and  if  the  amplitude  of  the  larger  event  is  more  than  three 
times  that  of  the  smaller  event.  This  is  also  true  for  measuring  a reliable 
Mg • This  agrees  well  with  previous  studies  of  FK  .uMB  vs.  maximum  likelihood 
filtering  (Lambert  and  Der,  1974 j. 

In  applying  the  method  it  can  be  assumed  that  the  events  in  question 
have  been  detected  on  the  short  period  instruments  and  therefore  their  great 
circle  azimuths  are  known  and  the  separation  has  to  be  performed  in  order  to 
measure  the  excitation  of  surface  waves.  It  is  also  possible  to  estimate 
the  azimuth  of  the  larger  event  by  one  of  f-k  spectra.  This  might  be  useful 
practice.  A-priori  information  about  multipathing  iron  the  epicenter  might 
also  be  used.  The  estimated  azimuths,  instead  of  great  circle  azimuths,  could 
conceivably  increase  the  effectiveness  of  the  separation  since  the  actual 
directions  of  energy  arrivals  could  be  included  in  the  model.  Thus  in 
applying  the  method  several  options  are  available  where  one  can  trade  compu- 
tational efficiency  and  speed  for  more  effective  separation.  In  cases  of 
complicated  multipathing  it  may  be  impossible  to  untangle  the  waves 

satisfactorily  and  M values  may  simply  have  to  be  made  from  the  f-k  spectra. 

s 
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Figure  2.  Polar  plots  of  energy  rejection  in  dB  for  iterative  beam  vs.  si  pic  bea 
Baja  California  event.  Only  the  A,  B,  and  C rings  of  USA  are  utilized. 
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Polar  plots  of  energy  rejection  in  dB  for  iterative  beam  vs.  simple  beam 
•vent.  All  LASA  LPZ  sensors  used. 
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Figure  5.  Polar  plots  of  energy  rejection  in  dB  for  iterative  beam 
Guatemala  event.  All  LASA  LPZ  sensors  used. 
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Figure  6.  Polar  plots  of  energy  rejection  in  uB  for  iterative  beam  vs.  simple  beam, 
Guatemala  event.  Only  the  A,  B,  and  C rings  of  LASA  are  utilized. 
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Figure  7d.  Location  of  256  second  overlapping  time  windows  used 
Baja  Caliiornia  event. 


Figure  7e.  Location  of  256  second  nonoverlapping  time  windows  used 
San  Juan  Province  Argentine  event. 
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Figure  5b.  Time  domain  analysis  of  the  Philippine  and  New  Britain  events  using 
the  reduced  LASA  array.  For  detailed  explanation  see  caption  under  Figure  8. 
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Figure  11a.  Time  domain  analysis  of  the  Philippine  and  New  Hebrides  events  using 
the  reduced  LASA  array.  For  detailed  explanation  see  caption  under  Figure  8. 
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Figure  lib.  Time  domain  analysis  of  the  Philippine  and  New  Hebrides  events  using 
the  reduced  LASA  array.  For  detailed  explanation  see  caption  under  Figure  8. 
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Figure  12a.  Time  domain  analysis  of  Baja  California  and  Guatemala  events  using 
the  full  LASA  array.  For  detailed  explanation  see  caption  under  Figure  8. 
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Figure  13a.  Time  domain  analysis  of  the  Baja  California  and  Guatemala  events  using 
the  reduced  LASA  array.  For  detailed  explanation  see  caption  under  Figure  d. 


Figure  13b.  Time  domain  analysis  of  the  Baja  California  and  Guatemala  events  using 
the  reduced  LASA  array.  For  detailed  explanation  see  caption  under  Figure  8. 
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Figure  14a.  Time  domain  analysis  of  the  Baja  California  and  Argentine  events  using 
the  full  LAS\  array.  For  detailed  explanation  see  caption  under  Figure  8. 
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Figure  15a.  Time  domain  analysis  of  the  Baja  California  and  Argentine  events  using 
the  reduced  LASA  array.  For  detailed  explanation  see  caption  under  Figure  8. 


ITERATION 


